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ABSTRACT 

We investigate the abihty of both LTE and Non-LTE models to fit the near 
UV band absolute flux distribution, /a (A), and individual spectral line profiles 
of three standard stars for which high quality spectrophotometry and high res- 
olution spectroscopy are available: The Sun (G2 V), Arcturus (K2 III), and 
Procyon (F5 IV- V). We investigate 1) the effect of the choice of atomic line list 
on the ability of NLTE models to fit the near UV band fx level, 2) the amount 
of a hypothesized continuous thermal absorption extinction source required to 
allow NLTE models to fit the observations, and 3) the semi-empirical tempera- 
ture structure, 7kin(logT5ooo), required to fit the observations with NLTE models 
and standard continuous near UV extinction. We find that all models that are 
computed with high quality atomic line lists predict too much flux in the near 
UV band for Arcturus, but fit the warmer stars well. The variance among inde- 
pendent measurements of the solar irradiance in the near UV is sufficiently large 
that we cannot definitely conclude that models predict too much near UV flux, 
in contrast to other recent results. We surmise that the inadequacy of current 
atmospheric models of K giants in the near UV band is best addressed by hy- 
pothesizing that there is still missing continuous thermal extinction, and that 
the missing near UV extinction becomes more important with decreasing effec- 
tive temperature for spectral classes later than early G, suggesting a molecular 
origin. 



Subject headings: stars: atmospheres, late-type individual (alpha Boo, alpha 
CMi) Sun: atmosphere 
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Introduction 



Short &: Lesterl (Il996l ) and IShort &: Lesterl (11994) prese n ted an assessment of the abihty 
of LTE modehng performed with the ATLAS9 code Kurucz (jl992a ) to fit the near UV band 
of the Sun and the Kl.5-2 III giant Arcturus (a Boo, HR5430, HD124897), respectively. 
The former paper included an assessment of the quality of various input line lists, and found 
that the line list that provided the best match to the detailed high resolution near UV solar 
spectrum allowed too much flux to escape in the near UV band, making the model too 
bright compared to absolute spectrophotometric measurements in the near UV. The latter 
paper made the case for a hitherto undiscovered source of continuum thermal extinction in 
the near UV band of early K giants on the grounds that LTE models and synthetic spectra 
computed with any line lists, including the most complete ones, predicted signific antly more 



f lux in the near UV band than was observed. In the same year as the latter paper, iBell et al. 



( 119941 ) also examined the ability of various line lists to provide a fit to the near UV band 
solar spectrum, and concluded that there may still be missing near UV extinction in models 
of the solar atmosphere. 



by Non-LTE (NLTE) modeling of Arc turus (jShort fc Hauschildtl 



The LTE studies of IShort fc Lesterl (119961) and IShort fc Lesten (|l994r) were improved 



20031 ) (hereafter Paper I) 



and the Sun (IShort fc Hauschildtl 120051 ) (hereafter Paper II) carri ed out with the PHOENIX 
multi-purpose atmospheric modeling and spectrum synthesis code dHauschildt fc BaronlEgggl ) . 
The results of Paper I and Paper II were that more realistic NLTE modeling only made 
the near UV band d iscrepancies found with LTE models by IShort fc Lesterl (119941 ) and 
Short fc Lesterl (119961 ) even worse as a result of NLTE over-ionization of Fe I partially lifting 
the "iron curtain" of near UV Fe I line extinction, thus making models even brighter in the 
near UV. That Fe I suffers NLTE over-ionization (or, more p r operly , LTE itnder-ionization!) 
in late type stars is not a new result (for, example, see lRuttenI ( 1l986l )). but Paper I and II pre- 
sented the first quantitative assessment of the effect on solar near UV fx levels with atomic 
models of Fe I and Fe II containing over 1000 atomic E levels combined, and accounting 
for non-local multi-species coupling due to over-lapping transitions (shared radiation fields) 
among dozens of relatively abundant chemical species. 

In this paper we investigate two ways of addressing the problem of atmospheric models 
of late-type stars bei ng too bright in the near UV band: 1) we revisit the suggestion of 
Short fc Lesterl (119941 ) that there may be missing continuum thermal extinction in the near 
UV, but on the basis of NLTE modeling rather than LTE modeling, and 2) we attempt to 
construct semi-empirical models of the atmospheric temperature structure Tkin(logr) that 
more closely match the observed near UV flux levels. We also attempt to map out the 
dependence of a possible source of missing continuum extinction by investigating its depen- 
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dence on stellar parameters among an (admittedly) small set of stars for which there is both 
high quality spectrophotometry that extends into the near UV band, and high quality high 
resolution near UV band spectra: namely, the Sun (G2 V), Arcturus (K2 III), and Procyon 
(F5 IV- V, a CMi, HR 2943, HD 61421). 

2. Near UV band /a(A) distributions 
2.1. Absolute spectrophotometry 

2.1.1. The Sun 



We have assembled four relatively recent independent measurements of the solar irradi- 
ance spectrum in the visible and near UV bands, which are presented in Fig. [1] The four data 
sets are most discrepant with each other in the near UV band, as described is greater detai l 
belo w, and, in the case of th e data sets of lNeckel fc Labd (119841 ). iBurlov-Vasiljev et al\ (119951 ) 
and iLockwood et al\ (Il992l ). this may be because the correction for air mass extinction is 
more difficult in the near UV than in the visible band. The four data g ets were re-sampled a t 
the sampling interval of the highest resolution data set, AA = 3.0 A ( Thuillier et a/.ll2003l ). 
and then each was convolv ed to the resolution of the rn ost coarsely sampled data-set, i? ~ 70 
at 3500 A (AA = 50) A ( Burlov-Vasiljev et a/.lll995[ ). For this aspect of our investigation, 
we are concerned with the fit of models to the overall f\ level on the intermediate-band scale, 
rather than to finer spectral features. 



Neckel fc Labd (Il984j ) presents a heavily cited construction of the irradiance based on 
measurements of the disk center intensity spectrum, = 1,A), combined with Fourier 
transform spectrometer (FTS) spectra and center-to-limb variation data obtained with the 
McMath solar telescope at Kitt Peak, and is the data set to which we compared PHOENIX 
models in Paper II. Their short wavelength data set spans the 3300 to 6300 A range in 
intervals of 10 A. 



Like iNeckel fc Labs! (Il984j ) , iBurlov-Vasiljev et all (119951 ) also constructed the irradiance 
in bands of A A = 50 A from /^(/i = 1, A) measurements made in the 3100 to 6850 A range 
with the 0.23 m solar telescope at Peak Terskol. They estimate the uncertainty in their 
measured fx values to be 2.5% in the near UV band. They note that their irradiance values 
are 8% higher in the 3100 to 4000 A range than those measured by iNeckel fc Labd (119841 ) 
(notwithstanding high frequency oscillations in the comparison at the wav elengths of the 
Fraun hofer lines) as seen in their Fig. 11, and they conclude that the values of lNeckel fc Labs 
(jl984j ) are too low. 
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Lockwood et all ( 19921 ) measured the irradiance in the 3295 to 8000 A range in A A = 4 



A bands with the 0.6 m Morgan telescope at Lowell Observatory using pinhole optics to 
reproduce a stellar-like image of the Sun, and calibrating their measurements against similar 
measurements of the standard star Vega. They estimate the internal error of their fx values 
to be ~ 2%. They also find that their values on wavelength scales wider than that of 
Fraunhofer lines in t he 3500 to 4500 A range are as much as 10% larger than those of 



Neckel fc Labd (119841 ). and this discrepancy can clearly be seen in their Fig. 6. 



Thuillier et al\ ( 120031 ) report space based measurements of the near UV and visible band 



irradiance made with the UV and VIS spectrometers of the SOLSPEC instrument that was 
part of the ATLAS payloads and the SOSP instrument on board the EURECA satellite. The 
UV spectra span the 2000 to 3600 A range in 4 A AA intervals, and in the 3400 to 8500 
A range in 10 A intervals. Their study includes comparisons between the SOLSPEC and 
SOSP instruments and other space-based measurements that were made at about the same 
time, and finds agreement among them at the level of less than ^ 5%. A special advantage 
of this data set in the near UV band compared to the others compiled h ere is that the 



correc tion for air mass extinction is obviated for space-based observ ations. iThuillier et aL 
( 19981 ) reports a comparison between t he visible band data described in IThuillier et all toO^ 
and the data of iNeckel &: Labd (Il984j ). and pa nel e) oi their F i g. 2 a nd their Table V also 
show that their values are ~ 4% larger than iNeckel fc Labd (Il984j ) in the 3500 to 4000 
A region. 



2.1.2. Arcturus and Procyon 



BurnashevI (119851 ) presents a compilation of f\ distributions of a large sample of stars, in- 



cluding 1) their own measurements made at the Crimean Astrophysi cal Observatory, 2) mea- 



surements from the "Sternberg spectrophotometric catalogue" (see iGlushneva et al\ (119841 ) 
and references therein) consisting of data taken from 1970 to 1984 with the 40 cm and 60 
cm telescopes of the Crimean Station of the Sternberg Astronomical Instit ute and the 50 cm 



telescope of the Fessenkov Astrophysical Institute, and 3) measurements of iKharitonov et al. 



(Il978l ) taken with the 50 cm Cassegrain telescope and spectrum scanner with a photo- 
multiplier tube taken at the Fessenkov Astrophysical Institute from 1968 to 1986. These 
data sets all generally cover the A range 3200 to 8000 A with AA = 25 A or 50 A, and 
have been r e-redu c ed to a common s pectr ophot ometric system (the "C hilean system") by 
BurnashevI (119851 ). IGlushneva et al\ (119841 ) and IKharitonov et al\ (119781 ) state that the in- 
ternal accuracy of their f\ measurements is ~ 3.5%. Arcturus and Procyon are among the 
stars for which data was reported in all three of these sources. 
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Additionally, for Procyon, which is the brightest of the three objects in the UV band, 
we were also able to obtain a mid-UV spectrum recorded with the International Ultraviolet 
Explorer (lUE) space based observatory in the "long wavelength" band (LWR09108 in the 
NEWSIPS lUE archive). This data set covers the A range from 1850 A to 3350 A, and 
provides a check on the ability of the models to fit the fx distribution at shorter wavelengths. 



2.2. High resolution spectroscopy 



Sun iKurucz et al\ ( 1l984l ) obtained a spectral flux atlas of the Sun with the Fourier Trans- 



form Spectrometer at the McMath/Pierce Solar Telescope at Kitt Peak National Observatory 
(KPNO) in the A range of 2960 to 13000 A. The spectral resolution, R, is ^ 800 000 in the 
near UV, but we have degraded the data set to R = 200 000 by convolution with a Gaussian 
kernel to reproduce a more typical "high resolution" stellar spectrum, and to aid comparison 
with synthetic spectra computed with a realistic sampling, AA. The spectrum is useful for 
the prese nt investigation as lo ng as R is sufficiently large for weak spectral lines to be fully 
resolved. IWallace et all ( 119981 ) independently obtained a full disk spectrum of the Sun with 
the same instrument and similar R value. We also include the latter in our comparison to 
assess the variance among data sets. We note that the latter spectru na is darker by about 
0.02 to 0.05/^ (where is the continuum flux level) than the atlas of iKurucz et al\ (119841 ) 
in the cores of some, but not all, strong features such as the CN band head at A3883 (see 
Fig. [5]) and the Ca II H and K lines, although the two spectra are in close agreement for 
weaker features, and in the pseudo- continuum. This non-systematic variance complicates 
the assessment of models based on high resolution spectra because it is precisely in the cores 
of strong features that the role of the continuous extinction manifests itself. 



Arcturus iHinkle et all ( l2000l ) have obtained a high resolution {R = 150 000) spectrum 



with the KPNO Coude feed telescope and echelle spectrograph. We note that the obs erved 
spectrum only covers wavelengths down to ~ 3730 A. Because IShort &: Lesterl (119961 ) con- 
cluded that a surprisingly large amount of extra near UV opacity was needed to fit the 
high resolution photographi c spectrum of iGriffinI (119681 ) . we also incorporate a comparison 
of the observed spe ctrum of iGriffiij (119681 ) (also of i? ~ 150 000 in the near UV) to that of 
( iHinkle et a/.ll2000l ). W e find t h at th e two observed spectra are in close agreement, except 
that in the spectrum of iGriffinI (119681) strong absorption features tend to be no t as d eep as 
those of the spectrum of IHinkle et all (120001 ). It should be noted that IGriffinI (119681 ) com- 
mented on the difficulty and subjectivity of establishing a continuum rectification (/^ level) 
in the blue and near UV band as a result of over-blanketing by spectral lines, particularly 
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in the vicinity of the Ca II H and K hnes (A3934 and A3968) as a result of the broad Van 
der Waals damping wings that affect f\ throughout the 3900 to 4000 A range. Given the 
difficulties of calibrating photographic data, and correcting it for background scattered light, 
as well as the u ncertainty in co n tinuu m placement, the discrepancy is not surprising. How- 
ever, given that IShort fc Lesten (119961 ) based their conclusions on their result that synthetic 
spectra have spectral absorption features that were too strong, it behooves us t o investigate 
the e xtent to which that conclusion may be weakened by comparison to the (IHinkle et al. 
2000l ) spectrum. 



Procyon iGriffinl (119791 ) have created a high resolution spectral atlas of Procyon in the A 



3150 to 7500 A range based on spectra obtained with the Coude spectrograph of the Mount 



Wilson 100-inch reflector. We note that, like the spectral atlas of iGriffinI (11968), the detector 



was a photographic emulsion. Based on the instrumental profile described in iMaeckle et al. 
(Il975h . we estimate that R 75 000 in the near UV band (FWHM ^ 0.05A). The same 
difficulty of continuum rectification in the near UV band that wa s discu ssed above for the 
Arcturus spectroscopy also applies to Procyon, as noted in iGriffinI (119791 ). 



3. Models 
3.1. Modeling 

We have computed model structures and synthetic spectra for the three program stars 
that explore a number of variations in the modeling. The details of the modeling are as 
described in Papers I and II. The same modeling procedure is used for all three objects, 
except that Arcturus is modeled with spherical geometry because of its lower surface gravity 
(log (7 < 2), whereas the Sun and Procyon are modeled with plane-parallel geometry. We note 
that our models are in radiative-convective equilibrium (RE), and, thus, do not exhibit chro- 
mospheric or coronal heating in their outer atmospheres, which is caused by non-radiative 
processes. Therefore, our models do not have inversions in the Tkin(logr) structure, and 
cannot reproduce near UV spectral features that are chromospheric in origin (eg. the Ca II 
H & K line cores). However, for A > 3000 A, the flux in the continuum, and in weak to 
moderately strong spectral lines, arises from layers of the atmosphere that are well below 
the depth of te mperature mini r num, log T^nm {Tmm ) , in semi-empirical chromosp heric models 



of the Sun (see iFontenla et al\ ( 120071 ) and references therein) and Arcturus (see lAyres et al. 
( 119941 ) and references therein). Therefore, as long as we avoid any attempt to interpret the fit 
to the cores of highly opaque lines where the flux arises from the outer atmosphere, then the 
^kin(logT) structure of our models should be relevant to the spectral features being studied. 
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Similarly, our models are static and horizontally homogeneous, and so do not account for 
the complex spatial and temporal structures on a variety of scales that have been observed 
in the solar atmosphere. With the exception of sunspots and the granulation pattern, these 
structures (plages, faculae, Ca II K bright points) exist in the chromospheric layers and are 
probably related to the non-radiative chromospheric heating mechanism, and, as mentioned 
above, our models do not account for the chromospheric Tkin(logr) structure in any case. 



3.2. Stellar parameters 



In Tabled] we present the values of the stellar atmospheric parameters used in our mod- 
eling. Our parameters for Arcturus and the Sun are those of Papers I and II, respectively, 
and justification for our choices can be found there. We re-iterate here that we have cho- 
sen to investig a te sol ar models that are based on the "tradi tional" abundan ces, [A/H], of 
Grevesse et al\ ( 1l992l ) rather than the newer abundances of lAsplundl (120001 ) derived from 
horizontally inhomogeneous models f or con si stency with the ho r izonta lly homogeneous so- 
lar atmospheric modeling of iKuruczl (Il992al ). lAufdenberg. et al\ (120051 ) recently carried out 
a simultaneous atmospheric model fit to the spectrophotometric f\ distributions and near 
IR band interferometry for Procyon, and included multi-component modeling with combina- 
tions of synthetic f\ distributions from ID atmospheric models computed with PHOENIX to 
approximately simulate the effect of horizontal inhomogeneity. They find best fit values for 
the stellar parameters of Teg in the range 6510 to 6515 K, \ogg equal to 3.95 ±0.02, metallic 
[^] equal to 0.0, a depth-independent value of the microturbulent velocity dispersion, Xt, 
equal to 2 km s~^, and a value of the mixing length used in the approximate treatment of 
convection, l/Hp, in the range 0.5 to 1.25, where Hp is the pressure scale height. (They also 
find a value of the radius, R, equal to 1.46 x 10^^ cm.) 



3.3. Spectrum synthesis 

In every case the modeling is self-consistent in that the atmospheric structure was re- 
converged before computing the synthetic fx distribution for that model. For comparison to 
high resolution spectroscopy, we also produced continuum rectified synthetic spectra by di- 
viding the synthetic f\ distribution by a corresponding synthetic pure continuum distribution 
computed without line extinction. The synthetic spectra were computed with a sampling, 
AA, ranging from 0.015 A at A = 3000 A to 0.04 A at A = 8000 A , corresponding to 
a spectral resolution, R, of ~ 250 000. This ensured that weak spectral lines were at least 
critically sampled. The spectra were then resampled and convolved with an appropriate 
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kernel to match the R value of the observed spectra to which they were being compared. 



3.4. Models 



Model LTE-big is computed in LTE usin g the "big" inpu t line list that includes the full 
set of spectral absorption lines described by iKuruca (Il992al ). including millions of atomic 
lines, largely of Fe-group ele ments, t h at wer e theoretically predicted, but have not been 
observed in the laboratory. iKuruczl (jl992al ) found that this big line list was necessary 
to provide sufficient line blanketing of the synthetic fx distribution in the near UV band 
of a solar atmospher ic model to provide a good match to the observed f\ distribution of 
Neckel fc Labd (119841 ). This is the line list that was used in the LTE modeling presented in 
Paper II. 

Model LTE^small is the same as LT15-big except that the input line li st is the "small" 
list that was used in the ATLAS m odels of Kurucz prior to the ex pansion of iKuruczl (jl992al ) 
described above (see, for example, Kurucz &: Peytremann ( 19751 )). This is the line list that 
was used in the NLTE modeling presented in Paper II, and is the only line list that is 
consistent with the NLTE models produced with PHOENIX. 

Model NLTE is the same as LTE-small except that the coupled NLTE statistical equi- 
librium rate equations and radiative transfer equation are solved self-consistently for the 
lowest one to three ionization stages of nineteen astrophysically important light metals up to 
and including the Fe-group elements. This is the same model as the one designated NLTEpe 
in Paper II, and the particular species included in the NLTE treatment, and the details of 
the atomic models used, are those described in Table 1 of that paper. 

Model NLTE-extra is the same as model NLTE except that the continuous mass ex- 
tinction co-efficient, has been enhanced ad hoc by a factor that differs for each star, as 
described in Section [XBl 

Model LTE-small-extra is the same as model LTE-small, with the addition of the same 
amount of extra extinction that was added to the NLTE-extra model. Model LTE-small- 
extra allows an assessment of NLTE effects in the presence of extra near UV extinction by 
comparison to the NLTE-extra model with all other variables the same, including the UV 
n\ sources. 

Model NLTE-cool is the same as model NLTE except that the kinetic temperature 
structure, Tkin(T5ooo)) has been adjusted, ad hoc, as described in Section [3T6l The Sun was 
the only object for which this variation was either feasible or necessary. 
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3.5. Ad hoc extinction enhancement 

For each of the program stars, we produce NLTE-extra models by enhancing, ad hoc, the 
value of the continuous mass extinction co-efficient, k^, in the wavelength range where the 
discrepancy with the observed fx distribution was most severe, which was 3000 - 4200 A m 
all cases. The k^(A) distribution was enhanced in the affected AA range by adding to /t^ the 
value of / X k^'^, where / x k^'""" is the continuum thermal absorption coefficient computed in 
LTE, and / is a A-independent ad hoc factor. The k^'""" expresses the continuum extinction 
caused by the thermal destruction of photons only, as opposed to scattering extinction; ie. 
k\ = K^'""" + cr^, where, cr^ is the continuum scattering extinction coefficient. As a result, 
this enhancement models the effect of an as yet unaccounted for source of extinction that is 
purely thermal in that it only destroys and thermally emits photons rather than scattering 
them. The value of / was chosen separately for each program star based on visual inspection 
of the comparison of the observed and synthetic NLTE fx distributions, and is reported in 
Table [21 We note that given the speculative nature of this study, it was not our intention 
to determine the precise value of / that minimizes a statistical measure of the goodness of 
the fx fit, but, rather, to investigate whether such an ad hoc enhancement can plausibly 
improve the fit of NLTE models to both the absolute fx distribution and the rectified high 
resolution spectrum. 

The models with enhanced were then re-converged before computing the synthetic 
fx distribution so that the models are internally consistent. We find that the thermal equi- 
librium Tkin(T) structure of both the Sun and Arcturus, for both LTE and NLTE models, 
is negligibly affected by k\ enhancement below A values of 4200A because the frac t ion o f 



the star's bolometric flux being transferred in the UV band is small. IShort fc Lesterl (119941 ) 
found the same result for Arcturus with LTE ATLAS9 modeling. 



3.6. Ad hoc Tkin adjustment 

We produce the NLTE-cool model by adjusting the Tkin(logr5ooo) structure of the NLTE 
model downward in value by an offset, ATkin(log rsooo), that has its maximum amplitude, 
ATkin(— oo), at the top of the atmospheric model, which effectively has logrsooo = —6, and 
decreases linearly in value with increasing standard continuum logarithmic optical depth, 
logTsooo, to near standard optical depth unity (logrsooo = 0), from around which depth 
the emergent visible and near IR band continuum flux (/^) arises. The temperature offset 
at each depth in the model is computed as 

ATkin( logrsooo) = mx (logrsooo - (-6)) + ATkin(-oo), 
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restricted to the range —6 < log rsooo < 0, where the slope, m, is given by ATkin(— oo) / (—6), 
and ATkin(— oo) < because we are always cooling the atmosphere rather than heating it. 
The logTsooo > depth range is where alterations to the Tkin(logT5ooo) structure would 
most clearly amount to changing the de facto T^g value of the model as determined by 
the measured continuum flux distribution, Therefore, to maintain some distinction be- 
tween perturbations of Tkin( log rsooo) for a model of given Tes and perturbations to the Tes 
value itself, we restrict Tkin alterations to the log rsooo < range. The coupled equation 
of state (EOS), hydrostatic equilibrium (HSE), and the NLTE statistical equihbrium rate 
equations, and radiative transfer equation, were reconverged for the resulting model, but 
with the temperature correction procedure disabled so that a self-consistent model with the 
new Tkin(logr) structure was computed, then a synthetic fx distribution for this model was 
computed. The value of ATkin(— cxo) was chosen by visual inspection of the fit of the result- 
ing synthetic fx distribution to the observed near UV band fx values. This procedure is 
effectively semi-empirical modeling of the atmospheric Tkin structure with the near UV band 
fx distribution as the sole diagnostic. For the Sun, ATfcj„(— oo) = —400 K. For Arcturus, 
the most extreme value of ATfcj„(— oo) that we were able to investigate was 1500 K, and we 
found that there was no value of ATfcj„(— oo) that would sufficiently reduce the computed 
near UV fx distribution without the Tkin (log r) structure becoming so extreme that the fit 
to the well observed visible band fx distribution was lost, and, for the most extreme values, 
the population re- convergence became numerically unstable. 



4. The Sun 

4.1. Spectrophotometric fx distribution 

Fig. [T] shows the comparison of all four observed solar fx distributions with with the 
computed fx distributions from models with various treatments of the thermodynamic equi- 
librium and K^. Table [3] presents, as a quantitative measure of the relative goodness of fit, 
the root mean square (RMS) of the deviation of the synthetic from the observed fx distribu- 
tion relative to the observed distribution, a, computed from interpolated pixel-wise relative 
deviations in the 3500 to 4000 A range, 

^' = Efm-rx)/rxr/N, 

where fl and are the synthetic and observed fx distributions, respectively, and is 
the number of elements of the common A sampling to which all fx distributions have been 
interpolated. The a values have been quoted only to the first distinguishing digit. 

As reported in Paper II, model LTE-big reproduces the observed fx relatively well, thus 
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confirming the results of tlie ATLAS modeling of lKuruczl (jl992a( ). Also, as reported in Paper 
II, the NLTE model significantly over-predicts the fx values in the near UV band. In Paper 
II we attributed the entire discrepancy between the fx values of the LTE-big and NLTE 
models to NLTE effects, particularly NLTE over-ionization of Fe I by the UV band radiation 
field. However, a factor tha t was neglected in Paper II is that their LTE model incorporated 
the "big" atomic line list of lKuruczl (Il992al ) {ie. their LTE model is the same as our LTE-big 
model), whereas their NLTE mo del only incorporates the "small" line list that pre-dated 
the expansion of iKuruca (Il992al ). Therefore, some of the difference between the LTE and 
NLTEpe fx distributions is caused by the much larger number of lines included in the LTE 
model compared to the NLTE model. 

The model atoms used in the calculation of the NLTE line extinction in PHOENIX must 
be consistent with the line list because any transitions not included in the NLTE statistical 
equilibrium rate equations are taken from the line list. Therefore, it is not currently possible 
to compute a NLTE model with the "big" line list with PHOENIX. To remove the choice 
of line list as a variable, we have computed our LTE-small model, which is more directly 
comparable to our NLTE model, and to the NLTEpe model of Paper II, than the LTE model 
of Paper II is. In addition to the requirement of consistency, there is some justification for 
investigating the fit of models computed with the "small" line list in that the "big" line list 
has been controversial. The "big" list includes millions of transitions that were theoretically 
predicted using an approximate treatment of the quantum mechanical atomic physics, and, 
although these lines improved the broad-band fit to the solar UV fx distribution, comparison 
to high resolution solar UV spectra indicate that many of the gf values and transition 
energies are inaccurat e, and that, mor eover, many of the lines in the "big" line list seem 
to not actually exist (IBell et al\ Il994l ). This very important point is also demonstrated 
again here in Fig. [2l In general, discrepancies between observed and synthetic spectra may 
be caused by errors in the A values in the l ine list for indiv idual lines that really exist. 
However, we note that the careful analysis by iBell et all (119941 ) lead them to conclude that 
the "big" list really does systematically contain many lines of detectable strength that are 
either unobserved, or are much weaker, in the observed solar near UV and blue spectrum, 
whereas the reverse is not true, and this is also consistent with our experience . Fig. [2] shows 



the co mparison of the observed rectif i ed hi gh resolution solar spectrum of IWallace et al. 



(1l998l ) (the spectrum of iKurucz et all (Il984j ) does not extend down into this A range) and 
synthetic spectra computed with the models LTE-big and LTE-small in the A3540 to 3570 
region. Examples of where the "big" line list leads to synthetic spectra with lines that are 
either too strong, or are altogether unobserved, are scattered through the near UV band, but 
this spectral range i llustrates a numb er of easily discernib l e exa. mples well. This point has 
been made by both iBell et al\ (119941 ) and IShort fc Lesterl (Il996l ). but is important enough 
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to be worth demonstrating here. On the this basis we proceed on the assumption that the 
"big" hne hst is less reahstic than the "small" line list. 

The LTE-small synthetic f\ distribution is also shown in Fig. [T], where it can be seen that 
it, like the NLTE synthetic fx distribution, also over-predicts the observed flux, although 
by not as much. Therefore, the near UV band discrepancy between LTE and NLTE models 
reported in Paper II is reduced when consistent line lists are used. We conclude that the 
effect of NLTE Fe I over-ionization, while still significant, is not as large as that reported in 
Paper II. 



In Paper II we conclu ded on the basis of comparing synthetic fx distributions to that 
of iNeckel fc Labd (119841 ) that NLTE solar atmospheric models over-predict the near UV 
band flux by about 10%. However, including other solar irradiance measurements in the 
comparison weakens this conclusion. Fig. [1] shows that the NLTE synthetic fx distribution 
lies at the upper end, but within the range of, the fx values span ned by the four observe d 
distributions, and is generally consistent with that measured by iLockwood et al\ ( 119921 ). 
except for a region around 3500 A. We also note that the LTE-big fx distribution, although 
consistent with the data of iNeckel fc Labd (119841 ). is found to lie at the bottom of the range 
spanned by the four data sets, suggesting that the LTE-big model actuall y under-predicts 
the near UV band flux, in contradiction to the conclusions drawn by both iKuruczl (jl992al ) 
and in Paper II. By contrast, the LTE-small model, being less bright in the UV band than 
the NLTE model, lies well within the range spanned by the four data sets. This suggests that 
the "big" hne list of iKuruczl (Il992al ) actually contains too much line opacity to reproduce 



the ob served near UV fx level, which is consistent with the conclusions drawn by iBell et al. 



( 119941 ) based in high resolution spectroscopy. As a result, given the range spanned by the four 
data sets, there is no observational basis for claiming that the, presumably more realistic, 
NLTE model fits better than the LTE-small model, and there is marginal evidence that it 
provide a worse fit. This speaks to the need for a more precise observational de termination 
of the UV band fx values. Interestingly, the conclusion of both IKuruczl (Il992al ) and Paper 
II that the LTE-big model provides the be st fit may be due to a n unfortunate conspiracy of 
errors: the observed near UV fx values of INeckel &: Labd (119841 ) are too low, if one were to 
assume that the other data sets are more accurate, and, at the same time, the "big" line list 
contains spurious extra near UV line opacity that produces an LTE synthetic fx distribution 
that is also too low, thus apparently providing a good "fit" to the data that under-estimate 
the true flux level. 
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As can be seen from Table [3], the NLTE-extra synthetic fx distribution, by design, provides 
closer agreement to the dataset of iNeckel &: Labj (119841 ). whereas the NLTE model with no 



additi onal extinction provides clo ser agreement to the data sets of iBurlov-Vasiljev et al. 



( I1995I ) and lLockwood et al\ ( 119921 ) . The data set of iThuillier et al\ (120031 ) does not discrim- 
inate clearly between the two models. This may be interpreted as marginal evidence for 
missing near UV band continuous thermal absorption opacity with an enhancement factor, 
/, of approximately 0.15, but, that conclusion depends on which of the observed f\ distri- 
butions one uses in the comparison. We note that the NLTE-extra fx distribution retains 
the good agreement with observed fx values at longer wavelengths in the visible band. This 
is because the fraction of Sun's bolometric flux passing though the 3000 — 4200 A region is 
small enough that the k\ enhancement does not affect the model Tkin(T) structure. 



For 3000 < A < 3600 A, the LTE-small-extra spectrum tracks the LTE-big spectrum 
closely. This is not surprising given that the value of / was tuned to compensate for the lack 
of line opacity in the "big" line list so as to restore the good fit to fx provided by the latter. 
However, we note that in the region of 3700 < A > 4200 A the LT15-small-extra spectrum 
falls significantly below that of the LTE-big model, indicating that the effect of extra k\ 
does not exactly mimic that of the additional line opacity in the "big" list. However, we 
note that we arbitrarily and somewhat crudely made / independent of A in the 3000 to 4200 
A region, so we should not expect the LTE-small-extra model to reproduce the spectrum of 
the LTE-big model throughout that entire range. 



The existence of additional near UV continuous extinction, if real, could, in principle, affect 
NLTE models by extinguishing the near UV radiation field that plays such an important 
role in, among ot her things, the well-known NLTE over-photo-ionization of Fe I (see, for 
example. Paper II. lRuttenI ( 1986 )). An assessment of the impact of / x on NLTE effects 
can be made by comparing the synthetic spectra computed with the LTE-small-extra and 
NLTI5-extra models in Fig. [H We note that, to a first approximation, the difference between 
the spectra of the NLTE and LTE-small models, and that between the NLTE-extra and LTE- 
small-extra models, is about the same throughout the 3000 to 4200 A range. We conclude 
that the magnitude of / x is not great enough to significantly alter the character of 
the NLTE effects on the solar fx distribution. This may not be surprising given the modest 
value of / (0.15). A more careful and detailed analysis of the impact of / x on the 
NLTE equilibrium of the solar atmosphere and on the detailed solar spectrum is beyond the 
scope of this investigation, but may be warranted if evidence mounts that the / x source 
really exists. 
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Fig. [3] shows the comparison between observed and synthetic fx distributions for the LTE- 
small, NLTE, and NLTE-cool models, and Table [3] also includes a values for the comparison 
of the NLTE-cool fx distribution to the four data sets. The best fit NLTE-cool model has 
ATkin{—oo) = —400 K, corresponding to a ATkin(logr5ooo) offset from the NLTE model with 
a slope of ~ 70 K in the — oo < logrsooo < range. The NLTE-cool r aodel is, by design, in 
better agreement than the NLTE model with data of iNeckel &: Labd ( 1l984j ). although none 
of these models is significantly discrepant with the observed fx values. However at longer A 
values in the visible band the NLTE-cool model significantly under-predicts the observed fx 
values, as is expected for a model with a Tkin(T) structure that is significantly cooler than 
standard solar atmospheric models. 

Fig. m shows the Tkin(log Tsooo) structures of the NLTE-cool and NLTE models, and, for 
co mparison, the Tkinflogrsnnn ) str uctures of the LTE s e mi-em pirical solar at mospheric models 
ofEo 



weger fc Muelleii (119741 ) and lGrevesse fc Sau vail (119991 ). The model of lHolweger fc Mueller 



(I1974J ) is based on fitting the synthetic high resolution fx spectrum computed with an 
LTE model with an ad hoc Tkin (logrsooo) structure to a variety of spectral features and 
to the center-to-lirnb vari ation of the intensity spectrum, />,(A,cos^), whereas the model of 
Grevesse fc Sauvall (119991 ) is based on an LTE analysis of the Fe I excitation equilibrium. 
Note that neither of these models show evidence of the well-known chromospheric Tkin(logr) 
inversion that sets in near logrsooo) = —4 because both investigators pointedly avoided 
spectral features that show evidence of chromospheric influence when choosing their fitting 
di agnostics. Interesting l y, the (logrsooo) structure of model NLTE-cool is closer to that 
of lHolweger fc Muellerl (119741) than that of the NLTE model, whereas the the NLTE model 
structure is closer to that of iGrevesse fc Sauvall (119991 ). 



4.2. High resolution spectroscopy 



Fig. [5] shows the c ompa rison between the observed rectified high resolution solar spec- 
trum of iKurucz et all (119841 ) and synthetic spectra computed with the LTE, NLTE, and 
NLTI3-extra models. We have chosen three regions that demonstrate how fx/ fx varies when 
Kx/i^x varies greatly over a narrow AA range due to a variety of different opacity sources, 
such as strongly damped spectral lines from a number of chemical species, and molecular 
band heads. The value of / is small enough that there is only s slight difference between 
the spectra of the NLTE, and NLTE-extra models, although there is marginal evidence that 
the latter gives rise to slightly more flux in very strong features, such as some strong atomic 
lines and the CN A3883 band-head. 



On balance, the fit to high resolution spectral features does not by itself provide com- 
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pelling evidence for the existence of missing near-UV band opacity. On the other hand, 
the addition of continuous opacity in an amount needed to improve the fit to the observed 
spectrophotometry does not significantly worsen the fit to the high resolution spectrum. 
Therefore, high resolution spectroscopy does not rule out such missing opacity. 



5. Arcturus 

5.1. Spectrophotometric f\ distribution 

Fig. [6] shows the comparison of the observed fx distributions to the synthetic spectra 
computed with models that explore the variation in the treatment of the thermodynamic 
equilibrium and the value of k^. Table H] presents the a values for each comparison computed 
as described in Section I^?T1 The synthetic surface fiux spectra. Fx, have been geometrically 
d iluted to the flux at Earth , fx, with an angular diameter for Arcturus, 0, of 21.0 mas 



( (IGriffin fc Lynas-Graylll999l )). which, as demonstrated in Paper I, provides a good match 
to the visible and near IR observed fx distribution. As reported in Paper I, both the LTE- 
big and NLTE models predict too much flux in the near UV band as compared to all four 
observed fx distributions, although the NLTE model is more discrepant with the data than 
the LTE-big model. As with the Sun, to remove the choice of atomic line list as a variable 
among the models, we have recomputed the LTE model of Paper I with the "small" line list 
that the NLTE model also necessarily uses, producing our LTE-small model for Arcturus. 
As expected, this model is less discrepant with the NLTE model, but, as a result, is even 
more discrepant with the observed values. As with the Sun, we note that in Paper I we have 
over-estimated the amount of NLTE UV band fx brightening by using inconsistent line lists 
in the LTE and NLTE modeling. 



The need for some adjustment to the model that reduces the UV band flux is even more 
obvious for Arcturus than for the Sun. We have computed a new NLTE model for Arcturus, 
NLTI5-extra, in which the value of / is unity, ie. the value of k^'""" is double that of a standard 
model. This model provides a close flt to the overall observed fx level, although the detailed 
distribution deviates locally from the observed one, most notably around 3400 A and 4100 A. 
NLT15-extra represents a large amount of missing thermal absorption opacity in the near UV 
band of early K giants. We note that the amount of extra continuous thermal extinction that 
we require to force NLTE m odels to agree wi t h obs erved near UV band spectrophotometry is 
consistent with the results of lShort fc Lesten (119941 ). who found on the basis of LTE modeling 
with the ATLAS9 code that k^'""" had to approximately double. 
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The LTE-small-extra spectrum consistently lies just below the NLTE-extra spectrum 
throughout the 3000 to 4200 A region, and, as with the Sun, shows a relation to the 
NLT15-extra spectrum that is similar to that of the LTE-small spectrum to the NLTE spec- 
trum. Again, without performing a detailed analysis of the effect of / x k^'""" on the NLTE 
model, we conclude to a first approximation that the extra extinction does not change the 
overall character of the NLTE effects on fx- We note that, in contrast to the solar case, the 
LTE-small-extra spectrum lies substantially below that of the LT&big model. For Arcturus, 
the value of / is so large that it more than compensates for the extra extinction provided 
by the extra lines in the "big" line list. As a result, it is much more clearly the case with 
Arcturus than with the Sun that the inability of standard models to reproduce the observed 
near UV flux cannot be addressed by choice of line list. 



Peterson et al. I (Il993l ) concluded that Arcturus has a non-solar abundance distribution 



with some a-process elements being enhanced by about a factor of two with respect to the 
overall scaled solar value. We comp uted NLTE atmospheric models and spectra with the 
abundances of lPeterson et al. I ( 1l993l ) and found that while there were detectable differences 
among the low resolution fx distributions computed with the scaled solar and non-solar mod- 
els, they were much smaller than the the discrepancy between the observed and computed 
spectra without / x k^'"*", and also much smaller than the difference between the flux spectra 
computed with LTE and equivalent NLTE models. We conclude that the near UV band fx 
discrepancy in Arcturus is not caused by inaccurate model abundances. 



As with the Sun, we attempted to construct a model, NLTE-cool, for Arcturus, that 
improves the fit to the observed near UV band fx values by ad hoc adjustment of the 
^kin(logT) structure. We were able to investigate models with ATfcj„(— oo) as large as -1500 
K, corresponding to a ATkin(logr5ooo) offset from the NLTE model with a slope of ~ 250 K 
in the — oo < logrsooo < range. The ATfej„(— oo) = —1500 K model, shown in Fig. [TJ has 
an upper boundary temperature at logrsooo = —6 of only ^ 500 K, and models with larger 
ATfcj„(— cx)) values are numerically unstable because the hydrostatic equilibrium solution 
has negative pressures at some depths. Therefore, the model of ATfcj„(— cx)) = —1500 K 
represents the most extreme Tkin(logr) structure that we were able to produce, and even 
it does not reduce the predicted fx values in the near UV to the observed ones, as can be 
seen in Fig. [71 We conclude that of the two modifications to the model proposed here, 
an additional source of can more effectively address the discrepancy with observations. 
Moreover, model NLTE-cool predicts significantly too little flux in the well observed visible 
band. We note that model NLTE-cool represents the extreme limit of so-called 1.5D modeling 
of horizontal thermal inhomogeneity in the atmosphere of Arcturus with components that 
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all have Tes ~ 4250 K. A semi-empirical component that is significantly cooler than an RE 
model is still unable to reproduce the observed near UV flux, even when this cool component 
has a weight of unity in the "1.5D model". Therefore, we conclude that 1.5D models with 
a mixture of hotter and cooler components that all have T^s ~ 4250 K cannot resolve the 
near UV discrepancy for Arcturus. 



We note that previous authors have also investigated non-radiative c ooling in th e outer 
atmosphere of Arcturus to address other observational discrepancies. iRyde et all (120021 ) 
observed H2O vapour lines in the mid-IR spectrum of Arcturus and found that the Tkin(logr) 
structure computed with a MARCS model in RE had to be depres s ed ad hoc by ~ 300 K 
in the logrsooo < —3.8 range. However, the model of iRyde et all (120021 ) has a boundary 
temperature of ~ 1800 K, 1300 K warmer than our NLTE-cool model, and comes into 
agreement with the RE model at logrsooo = — 3.8, at which d e pth o ur NLTE-cool model 
is still ~ 1000 K cooler than our RE model. [Wiedemann et al\ (119941 ) constructed a semi- 
empirical Tkin(logr) structure to account for the observed central flux of the CO Au = 1 
near-IR lines. They compared this model to a semi-empirical structure with a chromospheric 
^kin(logr) inversion fltted to chromospheric spectral diagnostics rather than to an RE model. 
The surprising depth of the CO Au = 1 lines seems to demand a signiflcant alteration 
to horizontally homogeneous semi-empirical models based on visible and near UV s pectral 
features, yet, even the CO based semi-empirical model of [Wiedemann et all (119941 ) has a 
boundary temperature of ~ 2400 K at a depth of logarithmic column mass density of —2.0 
(cgs units). This column mass density corresponds to a logrsooo value of ~ —5.9 in our 
model where the T^in va lue is less than 800 K, which, again, is over 1200 K cooler than the 



Wiedemann et all ( ll994j ) model. That our NLTE-cool model is so extreme compared to other 



semi-empirical models with cooled outer atmospheres based on other spectral diagnostics is 
further evidence that this model is unrealistic. 



5.2. High resolution spectroscopy 



Fig. E] shows the comparison b etween the observed high resolution spectra of both 



Hinkle et al.\ ( 120001 ) and lGriffinl (119681 ). and synthetic spectra computed with our LTE, NLTE, 



and NLTE-extra models for the same spectral regions that were shown for the flt to the Sun. 
As with the Sun, the strongly damped Fe I lines have slightly brighter wings as a result 
of NLTE Fe I over-ionization. However, again, the difference is almost negligible for many 
lines. By contrast with the Sun, the NLTE-extra model produces strong broad spectral 
features that are signiflcantly brighter, and in much better agreement with the observed 
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spectrum, than either the LTE or NLTE models. This agreement can be seen for strongly 
damped lines of both Fe I, Mg I, and the CN A38 83 band head. I t is si gnificant that we 
re-affirm with NLTE modeling the conclusion that IShort &: Lesterl (119961 ) drew with LTE 
modeling: namely, that the same enhancement of k\ that brings the near UV absolute f\ 
levels into agreement with the observed values, and also simultaneously brings the strong, 
broad features in the near UV high resolution spectrum into better agreement with data. 



6. Procyon 

6.1. Spectrophotometric f\ distribution 

In Fig. [H] we present the comparison between observed and modeled f\ distributions 
for Procyon. The synt hetic F\ spectra have bee n de-projected to f\ values at Earth using a 



value of of 5.404 mas lAufdenberg. et al\ (120051 ). The NLTE models agree more closely with 



the LTE models for Procyon than they do for the Sun or Arcturus, and, by contrast with 
the results for the Sun and Arcturus, both the LTE and NLTE models match the observed 
fx distribution well throughout the near to mid UV band, including the range observed by 
lUE. If anything, there is marginal evidence that the models predict too little flux in the near 
UV band. We find no evidence for the need for extra extinction, nor any other adjustment 
to the modeling, in the case of the hottest program star. We conclude that by Teg values 
of about that of the Sun, or slightly hotter, the near UV extinction is well reproduced by 
models, and do not present further discussion of Procyon here. 



7. Extinction sources 



Figs. [To] and [TT] show the dependence on the standard depth variable, logrsooo, of 
the most important thermal absorption and scattering contributions to k\ for the Sun and 
Arcturus, respectively, at a A value of 3500 A, near the center of the band for which we 
are proposing an additional continuous extinction source. We also include the log Tssoo scale 
among the abscissae so that an assessment of the value of near monochromatic optical 
depth unity can be made at the relevant wavelength. However, the two scales do not differ 
greatly for any star. For both tx scales, we present continuum extinction optical depth, due 
to both thermal absorption and scattering extinction. To control the size of extinction co- 
efficient output files, we have only tabulated and plotted the kx values for ten sample depth 
points spaced approximately equally at intervals of Alogr = 0.8. To simplify these figures, 
we have presented the total bound- free {h-f) thermal extinction due to all included metals. 
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K^'™ (labeled "M b-f"). In Figs. [121 [13] we show the logrsooo dependence of the contribution 
of the most important individual metallic species to k^'™ at the same A value for the two 
stars. Figs. [TH and [TJ] and Figs. [TH] and [T7] show the variation of the same quantities with A 
for a logrsooo value of 0.3, close to where the near UV continuum flux arises, with a sampling, 
AA, of 25 A. In all cases, we show the results of both LTE and NLTE calculations for those 
K\ sources that can be computed in NLTE by PHOENIX. We note that this is the first 
presentation of continuous extinction by source for a solar or late-type star, computed with 
thousands of transitions of all astrophysically important light metals and Fe-group elements 
treated in direct multi-level NLTE. 

As expected, for the Sun and Arcturus, both the LTE and NLTE PHOENIX calcula- 
tions confirm the long known result that the near UV band values are overwhelmingly 
dominated by H^ h — f extinction near continuum optical depth unity. For both the Sun and 
Arcturus, the H I 6 — / extinction just begins to rival that of H^ 6 — / at the deepest layer 
of the model (logrsooo > !)• For the Sun the total metal h — f extinction, when computed 
in LTE, and Thomson scattering both rival the H I 6 — / K3500 value at the highest layers 
of the model (logrsooo < while for Arcturus, Thomson scattering becomes dominant 
in the outer atmosphere for logrsooo < ~2. However, these layers are not expected to con- 
tribute significantly to the emergent surface flux distribution {Fx{t = 0)) in the near UV 
pseudo-continuum. The H~ b — f process is a well understood and accurately computable 
Kx source, and there seems to be little scope for belief that its value is erroneous by a large 
enough factor to explain the over-prediction of near UV flux by models. 

For those kx sources that are treatable in NLTE with PHOENIX, we note that, for 
the Sun, NLTE effects depress the value of the total metal b — f extinction at 3500 A for 
logrsooo < ~2 by as much as 0.5 dex by rsooo ~ 0, and increase the value of the H I 6 — / 
extinction at 3500 A for logrsooo < ~2 by as much as almost one dex by rsooo ~ 0. For 
Arcturus, the NLTE effects on both the total metal b — f and H 6 — / K3500 values are similar 
to that for the Sun, except that the NLTE depression of the H 6 — / extinction begins further 
out in the atmosphere (logrsooo < —3). All NLTE kx values approach the corresponding 
LTE value near the bottom of the atmosphere, as expected. 

Figs. [12] and [121 and [TH] and [TT] show that the NLTE depression of the total metal 
b — f extinction in the outer atmosphere is driven by a corresponding depression in the b — f 
extinction of many, but not all, of the most important individual metallic kx contributors. 
For the Sun, these are Si I, C I, and Fe I for rsooo > 1- For rsooo < 1; they are Fe I and Si I 
for an LTE calculation, whereas, for a NLTE calculation. Si I is still important, but Mg I 
rivals Fe I due to the large drop in Fe I extinction in the outer atmosphere as a result of the 
well-known NLTE Fe I over-ionization. For Arcturus, the important metallic contributors 
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are similar to those for the Sun, except that C I is not important at large depth, and both an 
LTE and a NLTE calculation show that Fe I and Mg I are important in the outer atmosphere, 
whereas Si I is not. Paper I and Paper II contain a detailed discussion of NLTE effects on 
the extinction of metallic elements for Arcturus and the Sun, respectively. However, for both 
stars, the effect on the total continuum extinction is almost negligible. Therefore, it does not 
seem likely that the over-prediction of near UV flux is due to the inaccurate or incomplete 
treatment of NLTE effects in continuum extinction sources (the effect of NLTE on line (b — b) 
extinction, especially that of Fe I, is very significant, as described at great length in earlier 
literature (see Paper II and references therein), but we are only investigating the role of 
continuum extinction here). However, as noted above, for none of our models does metallic 
b — f extinction dominate the near UV total value, so, although these NLTE results for 
continuum metal extinction may be interesting in their own right, it seems unlikely that 
inadequacies in their NLTE modeling will be able to account for the failure of NLTE models 
to correctly predict the observed near UV fx level in the Sun and Arcturus. 

Figs. [To] and [T^ show the logrsooo dependence at A = 3500 A , and the A dependence 
at logTsooo = 0.3, respectively, of the standard continuous extinction, k^, and the enhanced 
value, + / X K,'^^^ for the Sun, and Figs. [TT] and [15] show the same for Arcturus. Because 
the near UV extinction of both objects is dominated by the same sources, it does not seem 
likely that the missing kx source represented by / x k^'""" is directly related to any of the 
extinction sources shown in Figs. [TO] through [TT] le., it does not seem feasible that / x k^'""" 
represents an error in the calculation of any of the standard, well known sources because, 
presumably, an error of the magnitude required to fit the near UV fx level of Arcturus 
(/ ^ 1) would also have a significant impact on the Sun, whereas for the Sun / ^ 0.15. 



8. Conclusions 



One unanticipated conclusion of our i nvestigation is that th ere is a surprisingly large 
amount of scatter 10% if one includes the lNeckel fc Labsl (ll984j ) data set) in the measured 
absolute fx distribution for A < 4000A of the Sun. This undermines somewhat our ability to 
precisely test models of solar type stellar atmospheres, particularly with regard to the calcu- 
lation of the near UV extinction, kx- We conclude that a NLTE model based on the "small" 
high quality line list with standard extin ction is consistent with the observations, and pro- 
vides an especially good fit to the data of iLockwood et al. and lBurlov-Vasiljev et al. 



(Il995l ). The fainter near UV band data of iNeckel fc Labs! (119841 ) can be better fit if one 



adopts a modest enhancement, 0.15 x k^'""", in the 3000 — 4200A region. We are unable to 
compute NLTE models with the "big" line list, but we infer from a comparison of the LTE- 
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small and the NLTE synthetic spectra that s u ch a model would provide a closer match to 



the observed fx distribution of iNeckel &: Labd (Il984j ). thus obviating the need for additional 



continuum extinction that has been proposed previously on the basis of comparison to that 
particular data set. However, synthetic spectra computed with the "big" line list have been 
shown to be more discrepant to with the observed spectrum than those computed with the 
"small" list. The continuum-rectified high resolution solar spectrum is not sensitive to vari- 
ations in the value of of the magnitude necessary to reconcile the absolute fx distribution 
with observations. 

Interestingly, a NLTE model with standard extinction and an ad hoc "semi-empirical" 
^kin(logT5ooo) structure that is ~ 400 K cooler at a logrsooo of —6, our NL TE-cool model. 



i s also provides a better fit to the observed near UV fx distribution of iNeckel &: Labs 



( 1l984j ). and is in good agreement with a much more car e fully derived semi-empirical model 



^kin(logr5ooo) structure found by iHolweger fc Muellerl (Il974j ). However, the NLTE-cool 



model provide a worse fit to the well-observed visible band spectrum, and we do not be- 
lieve that adjustment of the Tkin(logr5ooo) structure of a horizontally homogeneous model is 
the most compelling resolution of the solar UV flux problem. 

The situation is much more clear for the red giant, Arcturus. All models, computed with 
any line list, grossly under-estimate the observed near UV band absolute fx distribution in 
the near UV. The standard must be approximately doubled to match fx in the 3000 — 4200 
A range. At the same time, the enhancement of also brightens the cores of strong 
spectra features and the wings of strongly damped spectra lines, thereby improving the fit 
of the synthetic spectrum to the observed continuum-rectified high resolution spectrum. We 
conclude that there is compelling evidence that there is still significant extinction missing 
from the models of early K star atmospheres, and that the missing extinction is probably 
due to a continuous process rather than line absorption. 

Given that the value of in the near UV in both the Sun and Arcturus is dominated by 
extinction sources that are readily calculated accurately, at least in LTE, namely, H^ b — f, 
and, to a lesser extent at extreme depths or heights. Hi b — f, and Thomson scattering, it 
seems unlikely that the discrepancy between observed and computed near UV fx levels is 
due to an inaccuracy in the extinction calculation of known sources, especially in Arcturus 
where the discrepancy to be accounted for is so large. Although NLTE effects are significant 
at height for both H I 6 — / and the total "metal" b — f extinction, the effect of NLTE upon 
the magnitude of is opposite for these two sources, causing NLTE deviations in the total 
extinction to cancel out significantly. Moreover, neither of these two sources dominates the 
near UV level around continuum optical depth unity. Therefore, it seems unlikely that 
the near UV fx problem is due to inadequacies in the treatment of NLTE alone, particularly 
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for Arcturus. However, we note that the computation of H~ h — f extinction only accounts 
for NLTE effects indirectly in the contribution of some, but not all, of the H and e~ sources 
in chemical equilibrium. 

By contrast, LTE and NLTE models of Procyon both reproduce the UV band f\ dis- 
tribution of Procyon down into the lUE LWR A range. If anything, there may be marginal 
evidence that models with standard extinction predict slightly too little flux in the near UV. 
Admittedly, these three stars provide a sparse set of data points for mapping out the depen- 
dence of any inadequacy in standard extinction calculations on stellar parameters. However, 
what is needed to improve the situation is a much larger set of standard stars for which there 
is both 1) high quality absolute spectrophotometry, /a(A), down to 3500, if not 3000, A, and 
high quality high resolution (sufficient to resolve spectral line profiles) continuum-rectified 
spectroscopy over a broad A range in the near UV, at least down to the limit for ground 
based astronomy {k, 3600 A). The sparseness of our stellar sample notwithstanding, com- 
parison of the ad hoc enhancement factor, /, for Procyon, the Sun, and Arcturus leads 
us to conclude that, if the near UV f\ discrepancy is due to an unaccounted for continuum 
extinction source, then that source becomes increasingly important with decreasing Tpff in 



the MK spectral class range from class mid-F to early K. We note that lShort fc Lester! (119941 ) 
did compare the measured f\ distribution of Arcturus with that of two other K1.5-K2 HI 
stars for which good quality data were available (but for which high resolution near UV 
spectroscopy was not available), and concluded that Arcturus is representative of stars of its 
spectral class in the near UV band. 



The dependence of / on T^q suggests a molecular origin, and IShort fc Lesten (119941 ) 
proposed the molecular photo-dissociation of metal hydride diatomic molecules that are 
relatively abundant in the atmospheres of late-type stars and that have lower dissociation 
energies than more common molecules, often corresponding to A values in the 3000 to 4000 
A range. Unlike photo-ionization, photo-dissociation leads to a broad k\ feature without 
a sharp edge because the upper, pre-dissociation, electronic state has a complex structure. 
However, this hypothesis is difficult to test quantitatively because values of the cross sections 
for relevant metal hydride molecules (MgH, SiH, FeH), either measured or calculated, are 
not reported in the molecular physics literature. 
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Table 1. Stellar parameters adopted. 



Object 


Teff (K) 


logg (cm s ^) 


m 


(km s-i) 


l/Hp 


Arcturus 


4300 


2.0 


-0.7 


2.0 


2.0 


Sun 


5777 


4.4377 


0.0 


1.0 


1.0 


Procyon 


6510 


3.95 


0.0 


2.0 


1.0 
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Table 2. Enhancement factor, /, of the ad hoc k\ enhancement used in the formula 

+ / X K'^ (see text). 



Object / 

Arcturus 1.00 

Sun 0.15 

Procyon 0.00 



Table 3. Sun: RMS (cr) values for goodness of fit of synthetic to observed f\ distributions. 



Data set 

Model NeckeP Burlov^ Lockwood^ Thuillier^ All data 



LTE-big 


0.048 


0.083 


0.082 


0.062 


0.071 


LTE-small 


0.040 


0.082 


0.084 


0.056 


0.068 


LTE-small-extra 


0.059 


0.118 


0.123 


0.093 


0.102 


NLTE 


0.074 


0.050 


0.039 


0.042 


0.053 


NLTE-extra 


0.027 


0.074 


0.076 


0.046 


0.060 


NLTE-cool 


0.048 


0.092 


0.093 


0.066 


0.077 



Neckel fc Labsl (|l98J) 



iBurlov-Vasiljev et all ( 119951 ) 



iLockwood et aLl fll992[ ) 



iThuillier et all (120031) 
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Table 4. Same as Table [31 but for Arcturus. 



Model 



Data set 

Burnashev^ Gluslineva^ Kharitonov"^ 



All data 



LTE-big 0.25 

LTE-small 0.31 

LTE-small-extra 0.20 

NLTE 0.46 

NLTE-extra 0.12 



0.32 
0.38 
0.20 
0.54 
0.15 



0.34 
0.41 
0.19 
0.56 
0.14 



0.30 
0.37 
0.20 
0.52 
0.14 



BurnashevI (119851 ) 



'^Glushneva et al. 


1984) 


^Kharitonov et al. 


(1978) 
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Fig. 1. — Sun: Comparison of observed and computed near UV band logarithmic flux 
distributions (log /a (log A)). Upper panel: Direct comparison. Lower panel: Difference 
bet\ veen observed ( f p and synthetic ifX) spectra (log.f°/.f? = log (f° — fXS) . Observed spec- 
tra: iNeckel fc Labd (119841): triangles: iBurlov-Vasiljev et al\ ( 119951 ): crosses: iLockwood et al. 



( I1992I ): squares; iThuillier et al\ (120031): Xs. Computed spectra: LTE: dashed lines; NLTE: 



solid lines. LT&big model (same as LTE synthetic spectrum of Paper II ( IShort &: Hauschildt 



20051 )): black dashed line; LTE-small model: dark gray dashed line; LTE-small-extra model: 



light gray dashed line; NLTE model: dark gray solid line; NLTE-extra model: light gray 
solid line. Note that, for simplicity, we only include logf^/fl for the NLTE and NLTE-extra 
models in the lower panel. 
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3560 3561 3562 3563 3564 3565 
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Fig. 2. — Sun: Comparison of observed spectrum of IWallace et al\ (1l998l ) (dotted line) with 
computed LTE high resolution flux spect ra (/a(A)) with various atomic line lists. (The 
observed spectrum of iKurucz et al\ (119841 ) does not extend down to this A range.) Model 
LTE-big: dark line; model LTE-small: light line. Note the relative quality of the model fits 
around Fe I A3549.7, the feature around A3550.8, around Ni I A3552.5, in the far red wing of 
the Fe I A3554.9 line, and the features (unidentified in the "big" line list) around A 3557.2, 
3557.7, 3559, 3561.5, and 3563.3, where the LTE-big model predicts spectral line absorption 
that is much stronger than that of either the observed spectrum, or the synthetic spectrum 
of the LTE-small model. The feature around A3550.8 is not identified in the "big" line list. 
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Fig. 3. — Same as Fig. [H except that the models explore the variation in Tkin(T5ooo) structure. 
LTE-small: darkest color; NLTE model: medium gray color; NLTE-cool model: lightest 
color. 
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Fig. 4. — Sun: atmospheric 7kin(log Tsooo) structure. Theoretical models: NLTE : dark line; 
NLTE^cool: lighter line. Semi-empirical models of iHolweger fc Muellen (119741 ): triangles; 
Grevesse fc Sauvall (1l999l ): squares. 
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Fig. 5. — Sun : Com parison of observed spectra of iKurucz et al\ ( 1l984j ) (dashed line) and 



Wallace et al\ (119981 ) (dotted line) with computed high resolution flux spectra (/a(A)). Model 



LTE: dark line; model NLTE: medium line, model NLTE-extra: light line. Note the relative 
quality of the model fits to the red wing of the Fe I A3749.5 line, Mg I A3829.4 and Fe I 
A3825.9 lines, and to the CN A3883 band head. 




3500 4000 4500 

log 10 ^ (A) 



Fig. 6. — Arcturus: Comparison of observed and computed near UV band logarithmic 
flux distributions (log /a (log A)). Upper panel: Direct co mparison. Lower panel: Dif- 
ference between f° ari d fl, log f° /f^. Observed spectra: iBurnashevI (Il985l ): diamonds; 



Glushneva et all (11984 ): triangles; iKharitonov et al\ (119781 ): Xs. Computed spectra: LTE: 



da shed lines; NLTE: solid lin es. LTE-big model (same as LTE synthetic spectrum of Paper 
II ( IShort fc Hauschildtll2005l )): black dashed line; LTE-small model: dark gray dashed line; 
LTE-small-extra model: light gray dashed line; NLTE model: dark gray solid line; NLTE- 
extra model: light gray solid line. Note that, for simplicity, we only include log f^/f I for the 
NLTE and NLTE-extra models in the lower panel. 
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Fig. 7. — Arcturus: Upper panel: Atmospheric Tkin(logr5ooo) structure. Theoretical models: 
NLTE: dark line; NLTE-cool: lighter line. Lower panel: Same as Fig. [6|, except that the 
models explore the variation in Tkin(T5ooo) structure. LTE-small: darkest color; NLTE model: 
medium gray color; NLTE-cool model: lightest color. 
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Fig. 8. — Arcturus: Com parison of ob s erved and computed high resolution flux spectra 
(/a (A)). Observed spectra: iHinkle et al\ ( l2000l ): dashed line. M odel LTE: dar k line; model 
NLTE: medium line, model NLTE-extra: light line, dashed line; iGriffiru ( 119681 ): dotted line. 
Note the relative quality of the model fits to the red wing of the Fe I A3749.5 line, Mg I 
A3829.4 and Fe I A3825.9 lines, and to the CN A3883 band head. 
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Fig. 9. — Procyon: Comparison of observed and comp uted near U V ba nd logarith- 
mic flux distributions (log fx (log A )). Observed spe c tra; o f iBurnashevI (Il985l ) : diamonds; 
Glushneva et all (119841 ): triangles; iKharitonov et al\ ( 119781 ) . lUE: crosses. Computed spec- 
tra: LTE-big model: darkest line; LTE-small model: dark gray; NLTE model: light gray; 
NLTI3-extra model: lightest color. 
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Fig. 10. — Sun: Depth dependence of the contribution to the total monochromatic continuous 
extinction, k^, at A = 3500 A from all significant sources. To limit the size of output 
files, only every tenth depth point in our 50 depth point model is plotted. Extinctions are 
computed in LTE (solid lines) and NLTE (dashed lines). (Note that only two of the major 
contributors, H I b-f and "metal" b-f, are treatable in both LTE and NLTE with PHOENIX 
(see text)). The dotted line shows the total extinction with the addition of the ad hoc extra 
continuum extinction, + / x k^'""" (see text). The vertical line indicates the Tsooo value at 
which the A dependence of is plotted in Figs. [H] and [161 
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Fig. 11. — Arcturus k^(t): Same as Fig. [TOl except for Arcturus. The vertical line indicates 
the r5ooo value at which the A dependence of Kc, A is plotted in Figs. [TSlandfTTl 
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Fig. 12. — Sun /t^'™(r): Same as Fig. [TOl except that the contribution to the total monochro- 
matic continuous "metal" 6-/ extinction, k^'™, from all significant metals is shown. 
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Fig. 13. — Arcturus K^'™(r): Same as Fig. [T2|, except for Arcturus. 



-42- 



. NLT Tot 



-1 - 



CP 

N 

E 
o 



C7> 
O 



- NLT M b- 



-2 - 



-3 - 



-4 




NLT H b- 



_L 



2500 3000 3500 4000 

X (A) 



LTE Tot 
H" b-f 




4500 



Fig. 14. — Sun K^(A): Same as Fig. [TOl except that the A dependence in the near-UV band 
of K^, at Tsooo = 0.3 from all significant sources is shown. The vertical hne indicates the A 
value at which the depth dependence of k\ is shown in Figs. [TU] and O 
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Fig. 15. — Arcturus k^(A): Same as Fig. [TH except for Arcturus. The vertical line indicates 
the A value at which the depth dependence of Kc,X is shown in Figs. [TT]and[T3l 
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Fig. 16. — Sun /t^'™(A): Same as Fig. [12], except that the A dependence in the near-UV band 
of the contribution to k^'™, at T5000 = 0.3 from all significant metals is shown. 
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Fig. 17. — Arcturus k^'™(A): Same as Fig. [161 except for Arcturus. 



